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ABSTRACT:  A  unified  framework  of  continuum  elasticity,  inelasticity,  damage 
mechanics,  and  fragmentation  in  solid  materials  is  presented.  A  free  energy  functional 
accounts  for  thennodynamics  of  elastic  deformation  and  damage,  with  kinetic  relations 
for  inelastic  rates  adapted  from  prior  literature.  Average  fragment  velocities  follow  from 
momentum  conservation,  average  fragment  sizes  are  deduced  from  an  energy  balance, 
and  size  and  velocity  distributions  are  constructed  following  entropy  maximization.  The 
model  is  further  developed  to  study  concrete  subjected  to  ballistic  loading  scenarios. 

INTRODUCTION:  Constitutive  models  of  dynamic  defonnation,  damage  evolution, 
and  fragmentation  of  solid  materials  are  needed  in  order  to  describe  complex  physical 
phenomena  occurring,  for  example,  in  solid  body  collisions  and  ballistic  impacts. 
Previous  research  efforts  have  often  relied  on  a  patchwork  of  separate  models  consisting 
of  an  equation-of-state  for  pressure,  plasticity  model  for  deviatoric  strength,  damage 
model  for  stress  reduction/elimination,  and  fragmentation  model  for  average  fragment 
size  [Johnson  et  ah,  1997;  Raftenberg,  1998;  Silling,  2001].  Such  patchwork 
implementations  may  lack  consistency  from  a  kinematic  and  thermodynamic  point-of- 
view.  The  aim  of  the  present  study  is  development  of  a  self-consistent  theory  accounting 
for  the  above  mechanisms,  specifically  amenable  to  brittle,  crushable  solids  such  as 
concrete,  mortar,  and  cinder  block.  Static  and  dynamic  properties  of  concrete  follow 
from  previous  investigations  [Holmquist  et  al.,  1993],  while  fragment  sizes  and 
distributions  are  predicted  via  an  approach  accounting  for  momentum  and  energy 
conservation  [Grady,  1982],  history-dependence  of  damage  accumulation  [Miller  et  al., 
1999],  and  maximum  disorder  principles  [Enghnan  et  al.,  1987;  Grady  &  Winfree,  2001]. 

PROCEDURES,  RESULTS  AND  DISCUSSION:  Fundamental  aspects  of  the  model 
framework  for  a  defonning  and  fragmenting  solid  are  presented  here.  The  kinematic 
description  begins  with  a  multiplicative  decomposition  of  the  deformation  gradient: 

FaA  =  dxa / dXA  =  FEaaFDaA  ,  (1) 

where  xa  and  XA  denote  spatial  and  reference  coordinates,  FEaa  is  the  recoverable 

elastic  deformation,  and  FDaA  is  the  irreversible  deformation  associated  with  defects  such 

as  micro-cracks,  voids,  dislocations,  or  shear  discontinuities  evolving  within  the  material. 
The  spatial  velocity  gradient  then  follows  directly  from  (1)  as 
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Lab  =  dxa/dxb  =  FEaaFE-lah  +  FEaaFDBaFD-wpFE-xph , 


(2) 


with  the  superposed  "U"  a  material  time  derivative.  Standard  balances  of  linear  and 
angular  momentum  apply.  The  Helmholtz  free  energy,  on  a  per  unit  intermediate 
configuration  volume  basis,  is  assumed  to  exhibit  the  following  general  dependencies: 

(3) 

where  the  elastic  strain  2Ee afi  =  FEaagahFEhp  -8ap ,  with  gab  the  spatial  metric,  0  the 

absolute  temperature,  and  a  dimensionless  set  of  internal  variables.  Stress-strain, 
temperature-entropy,  and  internal  dissipation  relations  are  then  deduced  as 


_  jE—\  j^Ea 


fj  =  -de\j/ , 


P°.LD’-aM<r?l)2:0.  (4) 


with  P'p  =  JEFE~laa^crab-if/Sab^FEhp.  Here,  aah  is  the  symmetric  Cauchy  stress,  JE  is 
the  elastic  Jacobian,  <//  is  the  free  energy  per  unit  current  volume,  fj  is  the  entropy  per 
unit  intermediate  volume,  and  LD/2a  =  Fd/Fd  lBa  is  the  purely  inelastic  velocity  gradient. 


The  framework  is  hereafter  specialized  to  brittle  crushable  solids,  with  FDaA  now 
associated  with  crack  opening  and  sliding,  as  well  as  pore  collapse  during  compression. 
Let  A  represent  the  cumulative  local  micro-cracked  area  per  unit  intermediate  volume 
V  ,  such  that  D  =  ^  oc  AVl/ 3 ,  and  restrict  0  <  D  <  1 .  Furthennore,  let  cp  =  ^  =  JD - 1 
represent  volume  reduction  upon  crushing.  The  specific  free  energy  density  then  is 
postulated  as 

v  =  \k\e‘i)2  +  r(n)+r(/>),(5) 

where  Kn  K2,  and  K,  depend  on  cp  and  D,  G  is  the  shear  modulus,  EEp  is  the  elastic 

strain  deviator,  F  accounts  for  surface  and  concentrated  energy  of  micro-cracks,  and  Y 
describes  the  specific  heat  content.  Deviatoric  plasticity  follows  from  the  potential  Q : 

K  =  [A1-d)  +  b(pIpSJ}+cGXIK)\^,  (6) 

where  sEh  is  the  spatial  deviatoric  inelastic  strain  rate,  3 A2  =  2 sEhsDba  for  stresses 
exceeding  the  elastic  limit,  6ah  is  the  deviatoric  stress,  -3  p  =  craa ,  and  A,  B,  C,  p0,  N, 
s0,  and  <j0  are  material  parameters  [Hohnquist  et  al,  1993].  In  the  present 
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implementation,  we  assume  small  elastic  shape  changes  and  material  isotropy  such  that 
the  inelastic  spin  may  be  neglected.  Porosity  and  damage  evolution  are  controlled  via 


|0  (p<  pc);(p>  Pl), 

\ap  ( PC<P<P /), 


D  =  kA([-7id(pI  p0)). 


(7) 


where  a  ,  k  ,  and  nD  are  positive  constants,  and  the  bracket  notation  (v)  =  x\/x  >  0  and 
lx)  =  OVx  <  0  .  Fragmentation  of  the  local  volume  element  is  assumed  to  occur  when  the 
damage  reaches  a  critical  magnitude,  D  =  DF .  At  this  instant,  a  spatial  energy  balance 
relating  intact  and  fragmented  configurations  is  written,  in  rate  form,  as 


d (U  +  Ec) / dt  =  d C*y,(u  +  e,  +  es )  +  Ec )/ dt ,  (8) 

where  U  =  JE~'  {y  +  ijO) ,  u  is  the  internal  energy  per  fragment,  e,  and  es  are  energy 
due  to  the  relative  linear  and  spin  momentum  per  fragment,  and  Ec  is  the  equivalent 

kinetic  energy  of  both  the  continuum  element  and  the  center-of-mass  of  the  fragment 
cloud.  Summation  is  implied  over  all  fragments  in  (8),  which  yields  the  average 
fragment  dimension  b.  Subscale  distributions  of  cumulative  fragment  number  per  unit 
volume  (AO  and  velocity  relative  to  the  center-of-mass  follow  from  statistical  mechanics 
principles  [Grady  &  Winfree,  2001],  i.e.  N  =  N0  exp(-A063)  where  N0  =Nf/V  is  the 

mean  number  of  fragments  per  unit  volume,  while  spatial-nonlocal  averaging  methods 
describe  fragments  larger  than  grid  sizes  in  numerical  implementations. 
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